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ABSTRACT 



Aims. The accretion-powered pulsar 4U 1626-67 experienced a new torque reversal at the beginning of 2008, after about 18 years of 
steadily spinning down. The main goal of the present work is to study this recent torque reversal that occurred in 2008 February. 
Methods. We present a spectral analysis of this source using two pointed observations performed by Suzaku in 2006 March and in 
2010 September. 

Results. We confirm with Suzaku the presence of a strong emission-line complex centered on 1 keV, with the strongest line being the 
hydrogen-like Ne Lya at 1.025(3) keV. We were able to resolve this complex with up to seven emission lines. A dramatic increase of 
the intensity of the Ne Lya line after the 2008 torque reversal occurred, with the equivalent width of this line reaching almost the same 
value measured by ASCA in 1993. We also report on the detection of a cyclotron line feature centered at ~37keV. In spite of the fact 
that an increase of the X-ray luminosity (0.5-100 keV) of a factor of ~2.8 occurred between these two observations, no significant 
change in the energy of the cyclotron line feature was observed. However, the intensity of the ~ 1 keV line complex increased by an 
overall factor of ~8. 

Conclusions. Our results favor a scenario in which the neutron star in 4U 1626-67 accretes material from a geometrically thin disk 
during both the spin-up and spin-down phases. 

Key words, accretion, accretion disks — binariesxlose — pulsars: individual (4U 1626-67 ) — stars: neutron — X-rays: binaries 
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1. Introduction 

4U 1626-67 is a Roche-lobe filling binary system, where mat- 
ter from the companion flows from the LI point into an ac- 
cretion disk (iRevnolds et al.L Il997h . This low mass X-ray bi- 
nary (LMXB) consists of a 7.66 s X-ray pulsar accreting 
from an extremely low mass companion ((0.04 Mo for / - 
18°, iLevine et al.L Il988h . Although orbital motion has never 
been detected in the X-ray data, pulsed optical emission re- 
processed on the surface of t he seco ndary revealed the 42 
min orbital period (Middleditc h et al.l ri981), confirmed by 
IChakrabartvi (Il998h . Most likely the binary system contains a 
hydrog en-depleted secondary to reac h such a short orbital pe- 
riod (Pa czvnski & Sienkiewiczl[l98lh . The faint optical counter- 
part {KZ TrA, y~17.5) has a strong UV excess and high optical 
pulse fraction (McClintock et al., 1977, 1980). 

The sign of the torque in 4U 1626-67 has r eversed in only 
two occasions since its discovery by Uhuru (iGiacconi et al.L 
Il972h . Monitoring of the source by the Burst and Transient 
Source Experiment (BATSE) on board the Compton Gamma 
Ray Observatory (CGRO) starting in April 1991, found the pul- 
sar spin ning down, i r nplying a changed s i gn in the accretion 
torque ( Wilson et aU. 119931: iBildsten et all |1994|) . It was esti- 
mated that the reversal occurred in mid- 1990. More recently, 
our daily monitoring of 4U 1626-67 with Fermi/GBM starting 
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in 2008 August found the pulsar again spinning-up rather than 
spinning-down. SwiftfBAT observations allowed us to cover the 
evolution of this second torque reversal. The transition took 
place at around MJD 54500 (2008 Feb 04) a nd lasted approx- 
imately 150 days (ICamero-Arranz et a n i2010h . 

A 48 mHz quasi-periodic oscillation (QPO) has been 
detected in t h e X-- ra y emission jS hinoda et al], 119901 : 
Kommer s et al 1 119981) . lOrlandinietaiT (1199 8) inferred a 
neutron star magnetic field in the range of (2.4-6. 3)x 10'^ G. 
To compute this m agnetic field range a source distance of 5-13 
kpc was assumed (IChakrabartvi Il998l and references therein). 
A ~37 keV absorption cyclotron feature was found in the 
0.1-200 keV BeppoSAX spectrum (' Orlandini etall Il998h . A 
m ore detailed s t udy o f this feature with Suzaku can be found 
in llwakiri et all (l2012h . The X-ray broad-band continuum was 
fitted with low-energy absorption, a blackbody, a power law 
and a high energy cutoff. 

lAngeUni et al ] (ll995h using the ASCA observatory discov- 
ered for the first time a strong emission-line complex centered 
on IkeV in the X-ray spectrum of 4U 1626-67. The strongest 
line was identified as hydrogen-like Ne Lya at 1.008 keV. The 
strength of the neon emission compared to the expected iron L 
complex implied a large neon overabundance. 

We present a spectral analysis of 4U 1626-67 with Suzaku, 
using two observations from 2006 March and 2010 September. 
This allows us to better characterize the recent torque reversal 
that occurred in 2008 February. 
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Table 1. Best Broad Band Fit Spectral Parameters" 





Observation 




a'' 

norm 


E'' 


E'' 

'^Fold 


'^CRF 

cr^ 


BBody 

nornigg 


Flux^ 


constant 
factor* 






400015010 
(Mar 2006) 


1.46(2) 


0.846(9) 
0.010(1) 


18.9(5) 


13(1) 


36 9+"" 

-^"•^ -0.8 

3.6(6) 


229 
127±13 


0.55(9) 


0.85(2) 
0.8(fixed) 


1.38(205) 


405044010 
(Sep 2010) 


0.6(2) 


0.981(3) 
0.026(1) 


21.6(7) 


15(2) 


37.8(8) 
4.3(6) 


0.511(8) 


1.58(4) 


1.21(4) 
0.8(5) 


1.20(206) 



"(PHABS (8 GAUSSIANS+GAUSSIAN+ BBODYRAD) +POW)*HIGHECUT ^=GABS)*CONST 
^ X 10-' cm"- 

Photon Index with normalization units in photons keV~^cm~^s~' at 1 keV. 
''keV 

The vakies for the cyclotron line depths are 14±4 and 18±4 for the Mar 2006 and Sep 2010 observations, respectively. 
^ with normsB-R^n /dTn, where dio is the distance in units of 10 kpc, varying from 0.5-1.3 (assuming the 5 -13 kpc distance range reported for this source). 

^ X lO"' erg cm"^ s"' (0.5-100 keV), based on the XIS 03 flux calibration; XIS 03 is a combined spectrum from XIS and XIS 3. 

^ Intercalibration factors (PIN/XIS 03 and GSO/XIS 03). The GSO/XIS 03 factor was fixed to the value obtained for the second observation since it could not be well 
constrained. 



2. Observations and Analysis 



The instruments on Suzaku (iMitsuda et al.L 
Imaging Spectrometer (XIS: iKovama et al 



2007h are the X-ray 



2007!) CCD detec- 



tor cov ering the 0.3-10 keV b and, and the Hard X-ray Detector 
(HXD; iTakahashi et"!!! l2007h comprised of the PIN diode de- 
tector (PIN) covering the 10-70 keV band and the gadolinium 
silicate crystal detector (GSO) covering the 60-600 keV band. 
The XIS has four separate detectors, XIS 0-3, with XIS 1 being 
a backside illuminated CCD and the others frontside. XIS 2 was 
lost due to a micrometeor hit in late 2006. 

The first Suzaku 4U 1626-67 observation occurred in 2006 
March (MJD 53797) and thus before the torque reversal in 2008. 
During this observation the spacecraft aimpoint was optimized 
for the XIS detectors, and Suzaku was run in a 1/8 window mode 
with 1 s integration time, an effective area of 128 x 1024 pixels 
and an effective exposure time of ~72.2ks. The second observa- 
tion in 2010 September (MJD 55444) was also centered on the 
XIS aimpoint but run in a 1/4 window mode, with 2 s integration 
time, an effective area of 256 x 1024 pixels and an effective ex- 
posure time of ~14.01 ks. For this observation only 3 XIS detec- 
tors were available, as we discussed in the previous paragraph. 
The Suzaku dat a were reduced with tools from the HEASOFT 
v6.11 package (lArnaudl \l99^ and CALDB release 20110804. 
Following the Suzaku ABC guid^l the first step of our analysis 
was to reprocess the data by running the FTOOL aepipeline for 
the XIS, the PIN and the GSO. The pipeline performs calibration 
as well as data screening. 

Due to thermal distorsion of the Suzaku spacecraft the opti- 
cal bench shows wobbling. The source's PSF image is not cen- 
tered on a fixed position on the XIS CCDs but appears bluiTed. 
As part of the standard data reduction this effect is partially cor- 
rected by adjusting the spacecraft attitude file. The attitude file 
can be further improved by applying the FTOOL aeattcor2: 
based on the algorithm described by Nowak et al. (2011) the 
mean event positions as a function of time are calculated and 
used to further adjust the attitude. While we obtain sharp images 
for our first Suzaku observation, the second one suffers from 
additional attitude variability that can be present in Suzaku ob- 
servations since a change in the attitude control system on 2009 
Dec 18. In this case the FTOOL cannot fully correct the character- 



istic "double images". According to the documentatiorQ, how- 
ever, and since conservative pileup exclusion regions were ap- 
plied (see below), the XIS spectra are not significantly affected. 

We used the FTOOL pileest to estimate the amount of 
pileup in the XIS images and disregarded circular regions with 
> 8% pileup fraction during the spectral extraction. This is a 
careful approach since due to the asymmetric shape of the PSF 
only few pixels outside of the exclusion region reach pileup frac- 
tions above 2 - 4%. For the first observation the pileup fraction 
remained below this limit in the entire source extraction region (a 
circle of ~ 70 radius centered on the source image). In the sec- 
ond observation pileup fractions of up to 15% were reached in 
the center of the source image so the extraction region consisted 
of a circle of ~ 130 radius minus a central circular exclusion 
region with a radius of ~ 35 . 

The primary tool for extracting data products (spectra, 
lightcurves) from XIS data is xselect, which is part of the gen- 
eral HEAsoft distribution. Events in the XIS detectors were in 
the 3x3 and 5x5 editing modes. With xselect we applied to 
the cleaned events filters with the good times intervals avoid- 
ing telemetry saturation. Source and background regions were 
created from the image. The filtered events were then used to 
create source and background spectra for the desired extraction 
regions. Products from the two editing modes were added. The 
XIS spectra were grouped close to the energy resolution of the 
detectors (see Nowak et al.L[201 lb . We created response matrices 
and effective area files with the xisrmfgen and xissimarfgen tools, 
respectively. The calibration of the XIS is uncertain around the 
SiK edge, thus events in the 1.8-2.0 keV range were excluded 
in the fitting. For more information about this issue we refer the 
reader to the "XIS Calibration around the Si K edge" web pagi 
(see also Iwa kiri et all 120121; iNowak et all 1201 U ISuchvet al 
201 1, for different treatments on this matter). 

For the PIN, we extracted spectra from the cleaned event files 
using hxdpinxbpi. This tool produces the dead time corrected 
PIN source spectrum as well as the PIN background spectrum. 
The non X-ray background events provided by the Suzaku team 
can be downloaded from the HEAS ARC archivfl The PIN spec- 
tra were grouped to have a signal-to-noise ratio > 10 in each 
energy bin, and we considered spectra between 14-70 keV. We 



http://heasarc.nasa.gov/docs/suzaku/analysis/abc 



' ftp://legacy.gsfc.nasa.gov/suzaku/doc/general/suzakumemo-2010- 
05.pdf 

^ http://heasarc.gsfc.nasa.gov/docs/suzaku/analysis/sical.html 
ftp://legacy.gsfc.nasa.gov/suzaku/data/background/pinnxb_ver2.0_tuned 
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1 10 100 1 10 100 

Energy [keV] Energy [keV] 



Fig. 1. The main panels show the Suzaku broad-band spectrum before (left panel) and after (right panel) the 2008 torque reversal. 
Blue points denote XIS data (0.5-lOkeV), PIN data points are in red (14-70keV) and GSO in purple (50-70keV and 50-100 keV, 
respectively). The residuals after fitting only the continuum model are shown in the middle panel, with no line emission and cyclotron 
line feature included. The bottom panel shows the residuals after fitting the line complex at ~1 keV, the ~6.5 keV iron line and the 
~37 keV cyclotron line (see Table[T]for a list of fitted parameters). A zoom of the region of the line complex at ~ 1 keV is also shown 
for clarity. 



point out that due to the changes in instrumental settings the PIN 
response matrices, from the CALDB database, are different de- 
pending on the epoch of the observation. 

The GSO spectra were created starting with the clean event 
files using the hxdgsoxbpi tool. The background events were 
downloaded from HEASARC0. Response files were taken from 
the CALDB database and the GSO correction ARF file was 
provided by the HXD instrument tearrQ The distributed GSO 
Non X-ray Background is grouped, so GSO spectra have to be 
grouped accordingly, using the spectral binning given by the file 
gsobgd64bins.dat available at HEASARCQ The GSO spectrum 
was restricted to the 50-70 keV range for the first observation, 
and to the 50-100 keV for the second observation. 



3. Results 

The ~0. 5-100 keV broad-band spectrum obtained using XIS (a 
combination of detectors and 3), PIN and GSO were fitted in 
XSPEC 12.7.0. We used a canonical model which allowed us to 



compare our spectral study with previous \yorks dAngelini et al 
'^avdo et al.'. 'I979t lOrlandini et all Il998t iKrauss et al 



1995 



2007 : Pa in et al., 2010, among others) and update the long-term 



^ ftp://legacy.gsfc.nasa.gOv/suzaku/data/background/gsonxb_ver2.5 

* http://www.astro.isas.ac.jp/suzaku/analysis/hxd/gsoarf/ 

^ ftp://legacy.gsfc.nasa.gov/docs/suzaku/analysis/gsobgd64bins.dat 



X-ray flux histor y of 4U 1626-67 relative to the flux mea- 
sured by HEAP 1 / chaki-abar&LeLaU[i99l|&Mss.eLaU 
ICamero-Arranz et al.. 2010). 

The continuum model fitted to both observations in- 
cludes low-energy absorption, a blackbody component, and a 
power law with a high-energy cutoff at ~20 keV (PHABS 
(BBODYRAD+POWLAW) HIGHECUT)*CONST. The con- 
stant factor (CONST) is used to account for intercalibration be- 
tween XIS, PIN, and GSO. The abundance was set to wilm 
(Wilms et al. 2000). Spectral fits for both observations are shown 
in Figure [T] A broad iron line near 6.4 keV (GAUSS) and a 
~37 ke V cyclotron feature (GABS) im proved the fit for both ob- 
servations (see also llwakiri et al.L 12012 ). Table[T]presents the fit- 
ted values for these broad-band parameters from the final model 
(which included the low-energy lines described below). 

Initial fitting results for the broad-band model for both obser- 
vations indicated evidence for a strong line complex at ~1 keV, 
shown in Figure [T] middle panels. Following Angelini et al. 
(1995), we added seven narrow Gaussian lines (listed in Table|2l) 
to our model. The line widths were fixed at the ASCA values 
(10 eV for all lines except the 1.362 keV line width was fixed 
at 57 eV). Table |2] lists the final fitted line parameters for both 
observations. Only five of the seven lines detected with ASCA 
were detected in the 2006 Suzaku observation, while six ASCA 
lines were detected in the 2010 Suzaku observation. For the tree 
ASCA lines not well detected in 2006 and the one in 2010 the 
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Energy [keV] Energy [keV] 

Fig. 2. The top panel shows the 8 Gaussian lines from the ~ 1 keV complex detected by Suzaku, before (left panel) and after (right 
panel) the 2008 torque reversal. Here the dotted lines in the model denote the components that were fixed. The residuals after fitting 
only the continuum model are shown in the middle panel, with no line emission included. The bottom panel shows the residuals 
after fitting the line complex with those 8 lines (see Table|2]for more information about these lines). 



Table 2. Fit Results on Line Emission' 



Expected Line Energy (keV) A Intensity* EW (eV) 

Observation (keV) Mar 2006 Sep 2010 Mar 2006 Sep 2010 Mar 2006 Sep 2010 Mar 2006 Sep 2010 

(dof=205) (dof=206) 



OHeo- 


0.568 


0.568(fixed) 


0.568(fixed) 






0.09(fixed) 


0.09(fixed) 


3.2(2) 


1.3(2) 


OLytt 


0.653 


0.661(8) 


0.66(2) 


59 


5.4'- 


3 0+' 


4(4) 


13(4) 


7(3) 


Ne Hea'' 


0.914 


0.930(8) 


0.929(8) 


105 


119 


1.3(3) 


7(1) 


8(1) 


14(1) 


Ne Lycf 


1.021 


1.025(3) 


1.031(2) 


878 


1699 


3.3(2) 


24(1) 


22.6(2) 


50.9(8) 


Ne Ylefi'' 


1.084 


1.084(fixed) 


1.13(3) 




39 


0.2(fixed) 


3(1) 


3.7(2) 


6(2) 


NeLy/? 


1.210 


1.210(fixed) 


1.25(1) 


6 


96 


0.2(fixed) 


2(1) 


1.8(2) 


9(4) 


Mg Heff 


1.34 


1.32(2) 


1.38(2) 


36 


31 


0.4(1) 


1.8+7 


4(1) 


5(1) 


Ne Recombo'^ 


1.362 


1.48(2) 


1.52(2) 


137 


154 


0.9(2) 


5(1) 


11.1(3) 


15(5) 


Fe K-fluor. line^ 6.4 


6.5(2) 


6.78(7) 


59 


30 


0.8-1 


1.3(5) 


42.2(3) 


33(1) 



" the widths of the low-energy lines were fixed to the ASCA values {10 eV for all but the 1.362 keV line with 57 eV) 
photons 10"^cm""s"^ 

'-' we note that, for the 2010 observation, the A.x~ for ^his line is not well computed due to the spectral model not being able to adequately fit this line (see Fig. [2] right panel). 
^ the abundances of the lines at 0.914 and 1.084 keV can also be consistent with a blend of several L-shell iron line emission (see Angelini et al. 1995). 

the Mg Ho- line emission occurs at 1 .47 1 keV (here not resolved, and thus not detected, by Suzaku), and therefore this Mg line may be contributing to this feature. 
^ for the iron line the widths were fitted (0.5(3) keV and 0.14(7) keV, respectively) 



intensities were fixed, since the eiTors could not be well con- 
strained. 

The broad feature at 1.362 keV could alternatively be ex- 
plained by a recombination edge. Fits to both observations using 
a model that replaced the 1.362 keV Gaussian line with a re- 
combination edge yielded similar 6x^ to our model containing 
the 1.362 keV line. The edge (and temperature) obtained for the 
2006 and 2010 observations are 1.45+0.05 keV (42+^° eV) and 

1.47+0.03 keV (45;';25 eV) (fitting with a Gaussian fine gives 
1.48(2) and 1.51(2)keV, see Table |2]l, respectively. A posteriori 
inspection of the residuals for both observations, suggested the 
presence of an extra Gaussian component at energy ~1.34keV. 
Therefore, up to 8 Gaussian components were included in the 



fit (see Fig. ^ Ax^ of 283 for 205 dof, and 248 for 206 dof, 
for the first and second observations, respectively). Final fit pa- 
rameters for this line at ~1.34keV, tentatively identified as the 
He-like Mg, are listed in Table |2] The A obtained by setting 
the flux of a particular line to zero in the best-fit model indi- 
cates the significance of each line. In addition, we fit the broad- 
band spectrum using different combinations of the XIS detectors 
yielding similar results. Other models accounting for emission 
from a hot diffuse gas, with line emissions coming from several 
elements, were also fit (meka and vmeka in XSPEC notation). 
We changed the elemental abundances while performing the fits 
with the meka and vmeka models, however they did not give rea- 
sonable fits. 
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A luminosity increase occurred between the two observa- 
tions, changing from La).5-\okeV) *2.4xl0^'^d^^^^, to 6xl0^'*d^^^^, 
erg s~'. However, the value of the center of energy of the cy- 
clotron line did not exhibit a significant change. This behavior 
resembles the case of the Be/X-ray binary system A 0535-1-26 
(Caballero et al., 2009), but is contrary to sou rces like HerX- 
1 and V0332-H53 (Pottschmi dt et al.L 120111 and references 
therein) or GX 304-1 (,Klochkov et al.Ll2012h r wh ich show pos- 
itive or negative L-E crf correlations (see also jBeckeretafl 
12012 *). On the other hand, the column density, A^h, dropped by 
a factor of ~3 while the black body temperature increased by a 
factor of ~2. The rest of the broadband parameters were consis- 
tent between the two observations. Absorbed fluxes in the 0.5- 
lOkeV, 2-lOkeV, 2-20keV and 2-50keV bands are 0.18(4), 



0.159(2), 0.295(3), 0.51(3) xlO"^ erg cm-^s"' and 0.49(3), 
0.38(9), 0.64(1), 1.3(9) xlO"'' erg cm-^g-', for the 2006 and 
2010 Suzaku observations, respectively. 

4. Discussion and Conclusions 

4.1. Spectral Continuum 

Our results show that two years after the torque reversal in 2008 
the luminosity has reached almost the same level as in 1977 (see 
Fig. |3]l, the black body temperature has increased to ~ 0.6 keV 
and the pulsar shows the same spin-up rate. The system's ap- 
pearance in this state can be explained in terms of accretion 
onto the neutron star from a geometrically thin Keplerian disk 
provided the mass-transf er rate is limited to M > lO'^gs ' 
(IChakrabartv et al.Lll997l) . 

The pulsar's transition from spin-up to spin-down in 1990 
was accompanied by significant changes of spectral parame- 
ters. The black body temperature decreased from ~ 0.6 keV to 
~ 0.3 keV an d the power law photon ind ex changed from ~1.5 
(Vauahan & Kitamoto, 199 7h to -0.6 (IVaughan & Kitamotol 
ll997;.Angelinietal.,. 1995.: lOwensetaLl 119971) . The transition 
period lasted about 150 days and the source luminosity during 
this time was only a factor of 2 smaller than its highest level 
observed in 1977 (see Fig. O. This indicates that the mass- 
transfer rate during the torque reversal in 1990 was in excess 
of 5 X 10'5gs-i. 

Our Suzaku observations show that the pulsar's behavior 
during the torque reversal in 2008 (from spin-down to spin- 
up phase) was different. While the transition period also lasted 
about 150 days and the black body temperature changed from 
~ 0.23 keV back to ~ 0.6 keV in 2010, the photon index re- 
mained almost unchanged. The value of the order of unity ob- 
served in 2006 and 2010 is the same as the photon index ob- 
served by several missions be fore the torque revers al during the 



spin-down phase, e.g. A5CA (lAngeUni et al ]ll995h. BeppoSAX 
(lOrlandini et al.Vl998'), Chandra and XMM ()Krauss_eLaLl|2007t 
[Schulz et al., 2001). The torque reversal in 2008 occurred at 
slightly lower luminosity level than in 1990 and was followed 
by a brightening of the X-ray source on a time scale of more 
than 2 years. 

As we have shown in a previous paper 
dCamero-Arranz et al.L 2010), the torque reversals in 4U 1626- 
67 are unlikely to be associated with alternation of the accretion 
flow geometry between the Keplerian geometrically thin 
disk (during the spin-up phase) and a spherical flow or a 
hot, geometrically thick disk (during the spin-down phase). 
The corresponding scenarios (see, e.g. Yi & Vishniac, 1999) 
encounte r major difficulties explai ning 150 day transition period 
(see also IWiiers & Pringlel 1 19991) and the low level of X-ray 



flux measured before and after the torque reversal in 2008. 
Furthermore, the mass-transfer rate during both of the observed 
torque reversals (> 5 x 10''' g s"') was above the critical mass- 
transfer rate at which the transition of the accre tion disk to the 
hot, geometrically thick state can be expected jNaravan & Yi 
1995). Finally, observations of the strong emission-line complex 
centered on 1 keV and the iron K-fluorescence line at 6.4 keV 
during both spin-up and spin-down phases indicate that the 
temperature of the material surrounding the neutron star does 
not exceed a few keV. Conversely, the temperature of the 
material in a thick disk is close to the adiabatic temperature, 
which for the parameters of 4U 1626-67 is in excess of lOkeV. 
Thus, our results favor a scenario in which the neutron star in 
4U 1626-67 accretes material from a geometrically thin disk 
during both the spin-up and spin-down phases. 

The magnetospheric radius of a neutron star accreting ma- 
terial from a geometrically thin disk can be evaluated as ^ 

K{f?/M{2GM„,y'^f'\ where fi = {1/2)B,RI^ is the dipole 
magnetic moment, B, the surface field and /?ns the radius of 
the neutron star, /c is a parameter accounting for the geome- 
try of the accretion flow, which in the case of disk accretion 
i s limited to k >~ 0.5 (iGhosh & Lambl [19791 \ma^ [19961 
Bozzoetal 1 [2009). For the star to be in the accretor state its 
magnetospheric radius should not exceed the corotation radius, 

'"coi = (CMns^'s /47r^) . Otherwise, the centrifugal barrier at the 
magnetospheric boundary w ould prevent t he acc retion flow from 
reaching the stellar surface (IShvartsmanill970h . This condition 
in the case of 4U 1626-67 is satisfied if the mass-transfer rate is 
M > Mmin, where 



M^in - 6 X 10'4gs-> Kl'^m-''^p-''ff4o 



(1) 



Here a:o.5 = k/0.5 is normalized to its minimum value, m is the 
mass and Pjj the spin period of the neutron star in units of 
1.4Mo and 7.7 s. The dipole magnetic moment of the neutron 
star, yU3o, is normalized to 10-"* Gem-'. According to our mea- 
surements of the energy of the cyclotron line feature with Suzaku 
(36.8^°^ keV) the surface field of the neutron star in 4U 1626- 
67 is ~ 3 X 10'^ G. This is c onsistent with previou s estimate 
of the field strength made by lOrlandini et alj (Il998h from ob- 
servations of 4U 1626-67 with BeppoSAX. Putting this value to 
Eq. ([TJ and taking into account that 4U 1626-67 is a persistent 
pulsar one finds that the mass-transfer rate has never decreased 
below 10'^ g s"' in the time period between 1977 and 2010. The 
rate of mass-transfer, therefore, may vary in the system up to an 
order of magnitude. It is remarkable that the observed X-ray flux 
(and, correspondingly, the mass accretion rate onto the surface 
of the neutron star) varies with similar amplitude (see Fig. [3]l. 
This indicates that either the value of k is indeed close to its 
lower limit, or the flux variations are governed by a different 
reason. In particular, one can also envisage a situation in which 
the magnetospheric radius of the neutron star during the low lu- 
minosity phase is close to the corotation radius (this implies k 
to be substantially larger than its minimum possible value) and 
the centrifugal barrier at the magnetospheric boundary reduces 
the mass transfer rate into the stellar magnetosphere (for discus- 
sion see Perna et al. 2006). The mass-transfer rate in the disk 
within this scenario varies insignificantly as the pulsar switches 
between the spin-up and spin-down phases. 

While the black body temperature during 2006-2010 in- 
creased from 0.23 keV to 0.51 keV, the source luminosity in- 
creased only by a factor of ~2. This implies that during the 
brightening phase the eff'ective area of the black-body source 
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has decreased by almost an order of magnitude. A high lumi- 
nosity and a relatively small size of this source suggest that it 
is located inside the magnetosphere in the vicinity of the neu- 
tron star The accretion process in this region is fully controlled 
by the stellar magnetic field and plasma flows within the ac- 
cretion channel. In this light, variations of the effective area of 
the black-body source can be associated with a restructuring of 
the accretion channel (i.e. change of both its cross-section and 
the height of th e shock at the base of the accretion column) . 
Studies (see e.g. 'Burnard e t all 1 1 983t iBasko & Sunvaevl 1 1 976t 
[Kulkar ni & Romanova, 200^ suggest that the parameters of the 
channel depend on the geometry and physical conditions in the 
accretion flow at the magnetospheric boundary, as well as on the 
mode by which the material enters the magnetic field of the neu- 
tron star Our observations give no evidence of any significant 
change in the accretion flow geometry during the brightening 
phase (see above). It, therefore, appears more plausible to asso- 
ciate the restructuring of the accretion channel with change of 
the mode by which the accreting material enters the magneto- 
sphere. 

It is widely adopted that the accreting material can en- 
ter the magnetosphere due to the Rayleigh-Taylor and Kelvin- 
Helmholtz instabilities or/and the magnetic reconnection. The 
Rayleigh-Taylor instability can be at work if the angular ve- 
locity of the accreting material at the boundary does not sig- 
nificantly differ from the angular velocity of the neutron star 
itself. The magnetic field configuration in the magnetopause 
in this case is close to poloidal. If this condition is not satis- 
fied the Rayleigh-Taylor instability can be suppressed by the 
magnetic field shear in the magnetopause (Bur nard et all Il983t 
Likh sanov & Pustil'nik, 1996). This indicates that the relative ve- 
locity between the magnetic field of the neutron star and the ac- 
creting material at the boundary can be one of the key param- 
eters determining the structure of the accretion channel. Even 
small fluctuations in the mass transfer rate may lead to signif- 
icant variations of this parameter if the ma gnetospheri c raditi s 
is close to the corotation radius (see e.g. iPerna et all l2006h . 
Another possibility is an amplification of the magnetic field at 
the i nner radius of the disk governed by the dynamo action (see 
e.g. iRappaport et al.L [2004). The disk in this case may switch 
into a sub-Keplerian state and decrease the angular velocity of 
the material at its inner radius. However, the modeling of the 
evolution of the channel, and its influence on the properties of 
the black body source within these scenarios, is very compli- 
cated and beyond the scope of the present paper 

4.2. Emission Lines 

The observed change of the equivalent width (EW) of the Ne LyP 
line at 1.021 keV during the brightening phase of the pulsar from 
22.6 to 50.9 eV, as well as the increase of its intensity a factor of 
~8 reflect a transition of the accretion picture to a new phase. 
This may indicate either changes of plasma para meters at the in- 
ner radius of the accretion disk dKii et al.L Il986i) or changes in 
the wind caused by interaction between the Ke plerian disk and 
the st ellar magnetic field at the corotation radius (iLovelace et al.L 
Il995h . As shown by Perna et al. (2006), the mass outflow rate 
can be higher during the spin-down phase and lead to a peri- 
odic circulation of material between the inner and outer radii of 
the accretion disk. This feedback results in a redistribution of 
mass in the accretion disk and can lead to variations of the mass- 
transfer rate in the system. If this is the cause of the observed 
luminosity increase after the torque reversal in 2008, the time 
scale of radial motion of the material in the disk is ~ 2 years. 
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Fig. 3. The X-ray flux history of 4U 1626-67 relative to the 
flux measured by HEAO 1, in the same energy band, from pre- 
vious works (Chakrabarty et al. (1997): circles; Orlandini et al. 
(1998): triangle; Krauss et al. (2007): stars; Camero-Arranz et al. 
(2011) unfilled squares) and the two Suzaku observations (blue 
diamonds) in the 0.7-60 keV band. 



This situation can be realized if the a parameter at the outer ra- 
dius of the disk is relatively small, 

. ~ 0.02 m^l- f ,r ( , (2) 

IWcms-i/ \2yrj \10»'cm/ ^ ' 

and, correspondingly, the density of material in the disk signifi- 
cantly exceeds its av erage value for X-ray pulsars (for discussion 
see e.g. lKingl,l2012l). Here Cj is the sound speed in the disk at the 
distance R. The slow turbulization of the material in the disk can 
also explain the small value of the st rength of the torque noise 
observed from 4U 1626-67 (Chakra barty et all, 199 7). 

In the present work, following lAngelini et alJ (Il995h . we 
have decomposed the strong emission-line complex centered on 
1 keV found in the Suzaku X-ray spectrum of 4U 1626-67 (see 
Fig. 121). In both Suzaku observations the strongest line is the 
hydrogen-like Ne Ly g line at ~ 1 .025 ke V, in agreement with the 
ASCA observations dAngelini et al The broad feature at 

~1.48 keV found in our observations is close to the expected 
energy of the H-like Ne recombination edge at 1.362 keV. The 
measured 57 eV broadening of the line at ~1.48keV suggests 
that we may be als o detecting the r ecomb ination continuum of 
Ne, as suggested bv lAngehni et al. I (fT99l for the ASCA obser- 
vation. We point out that the Mg Ha transition occurs at energy 
1 .47 1 keV (which has not been resolved by Suzaku), therefore 
this Mg line may be blended with the recombination contin- 
uum of Ne. Furthermore, with Suzaku an eighth line centered 
at ~1.34keV was detected, which we tentatively identify as the 
Mg Heo' (see TablelH. 

During the spin-down period, comparing the strength of the 
low-energy lines detected in the Mar 2006 Suzaku observation 
with previous measurements by ASCA, we see that the overall 
decrease is of the order of ~2.5 (see Table |2]). Be fore the 2008 
torque reversal. X-ray missions like BeppoSAX (Owens et alj, 
1997), Chandra (Schulz e t al., 2001; Krauss etal., 2007) and 
XMM dKrauss et all l2007h were able to resolve only up to 3 
or 4 lines of this complex. In general, the measurements of 
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the intensity of those lines are in agreement with our results. 
During the new spin-up period after the 2008 torque reversal, 
the flux of the lines increased in general by a factor of ~5 ex- 
cept the Ne Lya line at 1.021 keV that increased by a factor 
of ~8. This result is not surprising taking into account such a 
change in luminosity. Unfortunately, during a previous spin-up 
period there was only one measurement of the ~ 1 ke V complex 
by the Einstein observatory in 1979, and it was not resolved due 
to the lower data quality. The overall luminosity decrease asso- 
ciated with the 1990 torqu e reversal w as a factor of ~3 in the 
0.5-lOkeV ran) ^(AngelinietaD, 119951 and references therein). 
iKrauss et al.l (12007 ) found a decrease of the equivalent widths 
of the NeHea 0.914keV and NeLyo- 1.021keV lines during 
the spin-down period. The behavior of the EW for the O Hea 
0.568 keV and OLyo- 0.653 keV lines was more steady during 
the same period. With Suzaku we found a large increase of the 
EW of the Ne Lyo' 1.021 keV line after the 2008 torque reversal. 
The EW for this line reached almost the same value measured 
by ASCA in 1993. For the other lines the increase is less pro- 
nounced. 

The strong emission-line complex at ~ 1 ke V has been sim- 
ply explained as bein g due to a substantial Ne overabundance 
dAngeUni et al.Lll995l) . Angelini et al. (1995) raised the possibil- 
ity that 4U 1626-67 is embedded in a shock heated gas from the 
supernova remmant that created the pulsar, but as they point out, 
the absence of stron g lines from Si, S, a nd Fe is inconsistent with 
this idea. Instead, ISchulz et alJ (1200 ih proposed that this cool 
material originated in the accretion disk or the mass donor, so 
that its composition is then a clue to the nature of the companion. 
ISchulz et alJ (l200lh detected double-peaked emission lines and 
strong absorption edges of Ne and O in their Chandra HETGS 
spectrum, thus pointing to the companion in 4U 1626-67 being 
the chemically frac tionated core of a C -O or O-Ne-Mg crystal- 
lized white dwarf. ISchulz et al.l (1200 ih also concluded that the 
strong edges are du e to absorption in cool, metal-rich material 
local to the source. IKrauss et aP (1200 7) found that the double- 
peaked emission lines and the line ratios implied that they are 
formed in a high-density environment, suggesting that they arise 
somewhere in the accretion disk. Krauss et al. (2007) also noted 
that their value for the temperature of the line formation region 
(> 10* K) is characteristic of the highly ionized, optically thin 
outer layers of the accretion disk. We note that, in our work, sin- 
gle Gaussian line profiles provided an acceptable fit to the ob- 
served line shapes, and physical models like meka or vmeka did 
not yield good fits. 
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